Ground state fidelity and quantum phase transitions in free Fermi systems 
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We compute the fidelity between the ground states of general quadratic fermionic hamiltonians 
and analyze its connections with quantum phase transitions. Each of these systems is characterized 
by a L x L real matrix whose polar decomposition, into a non-negative A and a unitary T, contains 
all the relevant ground state (GS) information. The boundaries between different regions in the 
GS phase diagram are given by the points of, possibly asymptotic, singularity of A. This latter 
in turn implies a critical drop of the fidelity function. We present general results as well as their 
exemplification by a model of fermions on a totally connected graph. 



Introduction- Quantum matter at zero temperature is 
known to exist in different phases that cannot be contin- 
uously deformed into each other without giving rise to a 
singular behaviour of some physical quantity. This sin- 
gularity, which extends its influence over a finite range 
of temperature, is due to the competition between dif- 
ferent kinds of quantum fluctuations trying to order the 
ground state according to distinct and alternative corre- 
lation patterns. These phenomena, known as quantum 
phase transitions (QPTs) recently attracted a big deal of 
attention pj. In particular ideas drawn from Quantum 
Information Q] e.g., quantum entanglement, have been 
proven to be conceptually insightful as well as to provide 
effective tools to investigate QPTs |3|- 

In this paper we aim at further investigating the ap- 
proach to QPTs advocated in Ref. jj]. There it was sug- 
gested, inspired by the results of Ref. that the bound- 
aries between different quantum phases can be analyzed 
in terms of an extremly simple quantity: the overlap 
i.e., the scalar product, between the ground states (GSs) 
corresponding to slightly different values of the coupling 
constants. The critical points being characterized by a 
large, possibly discontinuous, drop of the overlap. More 
interestingly, in Ref. Q it has been shown that, for some 
models, the overlap contains information about finite-size 
scaling and universality classes. We shall use a quantity 
strictly related to the overlap, the fidelity 

We shall deal with general systems of free- fermions . 
This is a class of quantum models that are i) physically 
relevant ii) simple enough to allow for a rather compre- 
hensive analytical treatment and, at the same time iii) 
sufficiently rich to feature non-trivial GS phase diagrams. 
More specifically the class of fermionic systems we are 
going to consider is the one described by the following 
quadratic Hamiltonians 



H 



t 4a< 



L 

1/2 (clBijc] + h.c.) , (I) 



where: the q's (cj's) are the annihilation (creation) op- 
erators of L fermionic modes, A, B <E Ml(M.) are L x L 
real matrices, symmetric and anti-symmetric respectively 
i.e., A T = A, B T = —B. 



The Hamiltonian family and their ground states, 
can be directly parametrized by the generic real matrix Z 
and the defining A := (Z + Z T ) /2 and B := (Z T - Z) /2. 
From this perspective the space of coupling constants of 
Q is the L 2 -dimensional full matrix algebra Ml (R) even 
though most of the Z's will give raise to rather unphysical 
Hamiltonians e.g., highly non-local. In order to obtain 
physically meaningful models one might take, for exam- 
ple, A to be the adiacency matrix of some graph whose 
L vertices label the fermionic modes. 

The central quantity of our study is the GS fidelity || 



T{Z,Z) := |(* ; 



(2) 



where \^z) denotes the GS corresponding to the Hamil- 
tonian {T|l, In the case in which Z = Z + SZ we will 
write !F(Z,Z) — e~ s ( z K The whole idea of the present 
approach is that QPTs are characterized by a singular 
behaviour i.e., drop (enhancement), of T{Z) {S{Z)). 

The Ground State- Let us start by an explicit analyti- 
cal characterization of the ground state of the quadratic 
fermionic Hamiltonian Jy. To this purpose we use the 
fact that, for even number of particles and sites, the (un- 
normalized) ground state of Q is given by 01 



=cxp 1/2 



(3) 



where Cj|0) = 0, V« and G is a L x L anti-symmetric 
matrix, determined by gG + h = where g — f/2($ + 
h = f/2($ - *) G M L (R) and $(A - B){A + B) = 
A 2 $,A* = $(A - B), A = diag(A!,...,A L ) > [nj. 
The Aj's comprise the single-particle energy spectrum. 
From these equations, by recalling Z = A — B and by as- 
suming A (and hence Z) to be invertible one can immedi- 
ately obtain g = $(1 + A^ 1 Z)/2 and h = ^>(t-A^ 1 Z)/2 
where A$ := <f> _1 A"I>. From which, ifg is invertible, it fol- 
lows 



G 



K L Z - 1 
A^Z + 1 



T-l T 1 / 2 -T" 1 / 2 



T+l T 1 /' 2 + T- 1 / 2 



(4) 



where T := A^ 1 Z. The eigenvalue equation above can be 
rewritten as ZZ^ = A\, from which it follows that T is 
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the unitary part of the polar decomposition of Z = A$T 
[l2|. The operator g is not invertible iff —1 £ Sp(T), 
let be p its degeneracy. If p is even (odd) then it is easy 
to find out a canonical transformation, parity-conserving 
(parity-flipping) , mappin g T into a T' where all the — l's 
arc transformed into l's [l7| . The corresponding g' be- 
comes then invertible and one has a GS with the form 
©. This also implies that detT = (— l) p corresponds 
to the GS in the sector with parity (-l)^ 1 ^ = (-l) p . 
Since the parity is conserved by the Hamiltonian Q , its 
change in the GSs is always related to level-crossings i.e., 
first order quantum phase transitions, and these of course 
give rise to a discontinuous drop to zero of the fidelity. 
Now we are in the position of proving the following 

Proposition If the ground state of (JIJ is in the even- 
number of particles sector then it has the form 



I*; 



^=i[cos(^/2)|00)„,_, + sin(^/2)|ll) V! _ I/ ] 



where {e ±ie "}^ 



(5) 



Sp(T) and |0)„ (|1)„) denotes the 
vacuum (occupied) state of the new pairs of fermionic 
modes c±„ obtained from the Cj's by means of a suitable 
L x L unitary U. 

Proof. From Ref. 01 it follows that G can be brought 
by an unitary transformation U into the canonical block 
formG = U(0 L -2M®G D )U T , G D = i©JL x t v a\ v) . Here 
with o~ v ^ we denote a a y Pauli matrix acting on the two- 
dimensional span of two single-particle modes (conven- 
tionally labelled by v and —v) and t v G K — {0}. By 
allowing the first L — 2M (an even number) t v 's to be 
zero we can write G = UGdU t where the sum over v 
goes now from one to L/2. By redefining the fermionic 
modes via U and normalizing the state vector © \^z) '■= 
|*z)/||*z)|| one then gets |* z ) = ^MOO^ _„ + 
s v \\\) v - v ], with c := (l + tl)- 1 / 2 ,s v —Uil + tl)- 1 / 2 . 
From the block decomposition of G |l3| it follows that 
Sp(G) = {±^!/}^=i 5 on the other hand from Q one 
has Sp(G) = {£tan(0 M /2)}£ =1 where {e 1 ^} is the spec- 
trum of the unitary T, As stressed above we can assume 
T G SO L (E) and then T = e K where K G o L (R) i.e., it is 
skew-symmetric. One can now apply the decomposition 



of QJ to K i.e., K n = UKU T = 



( 



s can 



vanish) and obtain Sp(T) = {e ±%0 "}J =x . It follows that 
t v = tan(#„/2) and from this JSJ is obtained. □ 

To illustrate the formalism now we shall analyze the 
basic and elementary case of a pair of fermionic modes. 

Example 1. Two fermionic modes: A = ea z , B = 
iAa v . Then H = e{c\c x - c\c 2 ) + A(cJ4 + h- c -)> from 
which Z = ea z - iAo y , ZZ^ = (e 2 + A 2 )l + 2Aea x 
A = diag(|e — A|, |e + A|). By explicit computation 
of T = \Z\- l Z one finds that T is either a x or io*. 
Of course this model is trivially soluble with eigenvec- 
tors {|10>, |01), |00) + |11), 1 00) - |11)} and eigenvalues 
{e, — e, A, —A}. The zero-eigenvalue lines e = ±A splits 
the e — A plane in four regions, in each of which the 
ground-state is given by one of the states above. By 



considering instead the two fermion system with Z = 
u x + pt — iAo v (see the complete graph described be- 



low) one finds that T = a x (T = e 



= tan" 1 A/fi) 



for /i 2 + A 2 < 1 (pi 2 + A 2 > 1). The unit circle dividing 
the two regions is the set of points where Z is singular. 

Example 2. Two fermionic modes: A = el, B = iAa y . 
Then H — e(c\c\ + ejea) + A(cJcJ + h.c), from which 
Z = el - iAo-y, ZZ^ = (e 2 + A 2 )l, $ = 1, A = 
Ve 2 + A 2 1, T = exp(-i0o- v ) where := tan -1 (A/e). 
Therefore G= {e~ i6aV -t){e- l9ljV +1)" 1 = -i tan(§ a y ). 
At variance with the former case now we have zero single- 
particle eigenvalues only in the trivial case e = A = and 
T{&) is a smooth family of commuting matrices. This 
example can be readily extended to L = 2M fermionic 
modes by considering A = ©^^e^l,,, B = i 
A v a v v) . Then one finds T = ®^ =1 exp 

tan-^A^), 

The fidelity- We give now an explicit evaluation of the 
of the fidelity I n order to do that we first recall the 
following result from Ref. [l4| 



M 



det(l + GtG) 1 / 2 



det(l + GtG) 1 ^ det(l + G^G) 1 ^ 



(6) 



Proposition 1 



2V2 



det(T + f)| 



1/2 



L/2 

nicos(e„/2)|. (7) 



where the second equality holds for T 1 T G SOl(R) 
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and Sp(T- 1 T) = {e ±m »}%\. When det(T- x T) 
one has T(Z, Z) = 0. 

Proof. To prove the first equality it is a straightforward 
calculation from Eqs. JBJ, (@} and | det T| = |detT| = 
1. Of course if -1 € Sp^" 1 ! 1 ) then T{Z,Z) = 
but this happens ncessarily when det(T^ 1 T) = —1. If 
T- X f e SO L {R) one has Sp(T- 1 f) = {e ie "}£ =1 = 
jgiie^jV 2 ^ ^g rg |. e q ua my from unitarity, second one 
from reality and speciality) thus one finds J-(Z, Z) = 

3^072 = UiLl I c<»(e v /2)|. 



□ 



When [T, f ] = one has 6^ = 0, (fj, = 1, . . . , L - 
2M), Q v = ±(0 U -0 V ), (v=l,..., M); therefore one gets 
the following result extending the analogous one obtained 
for the XY model |j] : suppose now that Z and Z are such 

that [T, f] = 0, then T{Z,Z) = \[f =x cos(^f^). For 
example the Hamiltonian of Sr. J, for all possible choices 
of {e}„ and {A^}^, gives rise to commuting T's . 

To summarize: the conceptual path to find the fidelity 
between two different ground states of JJJ is described 
by the following chain of maps Z,Z ^ M L (R) — ► T, f G 
Ox,(R) — ► Sp(T- x f ) C S 1 . The first arrow denotes the 
passage from the matrices Z, Z to the orthogonal matri- 
ces T and T i.e., their polar parts, the second denotes 
the diagonalization of T _1 T. While the positive part of 
Z i.e., \Z\, contains the single-particle energy spectrum, 
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in its polar part T are encoded all the data to defining 
the many-body ground state JSJ. In the points where Z 
becomes singular i.e., zeroes in the single-particle spec- 
trum, the polar part is not uniquely defined: by crossing 
these points one can have a non smooth change of the 
T, which in turn results in a non smooth change of the 
ground-state structure. If singularity of Z is achieved 
for finite (infinite) size L one has a first (higher) order 
QPT. This kind of singular behaviour is reflected in Eq. 
(|7|) which clearly shows that an abrupt decrease of the 
fidelity, moving from Z to the neighbouring Z ', can be 
caused by the appearence of an eigenvalue ofT~ 1 T with 
a 'large' negative part; the closer this latter gets to — 1 
the smaller the fidelity. 

An instance of this phenomenon is of course provided 
by the discontinuous "first order" QPTs of Example 1 
in that det(T" 1 f') = -1 when T and f are in different 
phases. A less trivial one can be obtained from the multi 
mode case in Ex. 2. Suppose that, for some one 
finds e.g., in the thermodynamical limit, that A„ t— > + 
and then 6 Vo = lim A „ ^o cos" 1 ^/^ + 4J 1 / 2 ) = 
(7r/2)(l — ej/ /|e„ |). If e„ can be driven (in the same 
limit) throughout zero one observes a jump of 9 Uo with 
amplitude 7r; this in turn implies —1 is in the spectrum 
of T- X f [T = T(e„ = 0+) and f = f{e Uo = 0")]. This 
is the mechanism responsible for the fidelity drop in the 
XY model observed in the paramagnetic-ferromagnetic 
QPT (l5|. The important aspect is that the drop of 
corresponds here at e 2 + A 2 Q = i.e., the vanishing 
of one of the single-particle energies, or, equivalently, a 
gaplessness in the many-body energy spectrum. This kind 
of correspondence between the fidelity behaviour and ze- 
roes in single particle spectrum can be further unveiled 
by using the following " perturbative" result 

Proposition 2 If T(A) = e K( - x \K <E o L (R) and Z = 
Z{X),Z = Z(X + SX) one has T{Z, Z) = e~ s ^ where 

S(Z) = -—Ti(K'SX) 2 + 0(5X 3 ) (8) 
16 

where K' := dK/d\. 

Proof. From Eq. Q), using AetA = exp(TrlogA), 
one has T{Z,Z) = exp{(l/2)Tr log[(l + T- 1 f)/2]}, by 
using the Taylor expansion of the logarithm one obtains 

e SW , S(X) = En=o L n where L n = Iz^UrfT-^T)- 
Now one can write T(X) = e K(x \K € o L (R) then ST = 
f -T = K'TdX + (1/2)[K"T + {K') 2 T]8X 2 + 0(SX 3 ) 
(prime denotes derivation). From unitarity one has that 
K', K" are traceless, hence L x w 1 / >8Tr(K '5X) 2 and L 2 w 
— 1/16(K'8X) 2 . From here using the above expansion one 
obtains Eq. ©. □ 
To illustrate this result let us consider the system in 
Example 2 i.e., K(X) — i 0uW a f u ) with e„ and A„ 
functions of the parameter A. In this case, by taking 
the derivative of Q v = t&n^ 1 (A lJ /e v ) one gets Tr(K') 2 = 
-2Y, v {0' y ) 2 and hence T{Z,Z) = e -s^s 2 +o(s\ :i ) ^ whcrc 

^2 - l/8L(^z) 2 with D» z = e^Al - e~ 2 A v e' v . 



This equation suggests that if some of the \Z\ eigenvalues 
ye^+Aj is zero, or asymptically vanishing for L — > oo, 
then the function S2 should have a sharp increase and 
accordingly the fidelity a sharp decrease. 

Complete graph- In the remaining of the paper we will 
provide a GS fidelity study for the fermionic system 
with an underlying topology of a complete graph. A de- 
tailed analysis of this diagram will be reported in a sep- 
arated publication |17|. More precisely, now we consider 
the Hamiltonian defined by the following matrix data 
= 1 + (//- 1)% , B(7),j = 7 si g n (i - i) where 
(i, j = 1, . . . , L). We analyze now some simple cases in 
the (/x, 7) plane. 

i) For 7 = the corresponding number-conserving 
single-particle Hamiltonian is readily diagonalized: 
SpA(fi) = {L + fi — l f fi—l}, where the first (second) 
eigenvalue has degeneracy one {L — 1). We see that 
/j = 1 corresponds to the vanishing of the lowest eigen- 
value (for any L). ii) For (/i = 0, 7 = 1) the matrices 
Z and Z^ = Z T become lower and upper triangular re- 
spectively. By explict computation one finds (ZZ^)ij — 
4min(L — i, L — j), which has the last column and row 
identically vanishing. Accordingly € Sp|Z(0,l)| V L. 
iii) Changing the sign of 7 simply corresponds to trans- 
forming Z into Z*. Since Sp(ZZ^) = Sp(Z^Z), the single 
particle spectrum is not affected by the transformation. 
Furthermore, since T(- 7 ) = T%) and detT = detT*, 
the overlap behaviour is symmetric with respect to the 
7 = axis. 

The GS fidelity in the ^-7 plane is reported in Fig.^f 01 ' 
L = 400. By using T min = vwa[T{Z,Z &ll ),T{Z,Zs 1 )\, 
where Zg^ = Z{ji + S^i,j) and Zs-y = Z(n,j + 5~f), i.e., 
by plotting the minimum of the fidelity with respect to 
variations in both directions, one gets a clear diagram 
of the 'degree of orthogonality' of neighbouring ground 
states in the parameter space. The drops in the GS over- 
lap identify all the QPTs of the system. In general, the 
actual value of the drop at a given point in the /1-7 plane 
depends both on the chosen variations <5/i, £7, and on 
the system size L. For even values of L and in the ther- 
modynamic limit, the phase diagram is reported in the 
left panel of Fig. [21 The boundary of the region /j, < 1, 
\j\ < 1 corresponds to a first order QPT. Indeed, while 
the ground state inside this region contains an odd num- 
ber of fermions (detT = —1), it becomes even outside 
(detT =1)- Such a QPT is simply interpreted in terms 
of a level crossing in the energy spectrum. The phase 
diagram of this first order QPT can be clearly calculated 
also for finite values of L. In the right panel of Fig. [21 
the cases L = 2,4, 6,00 are shown. The line [i = 1, for 
I7I > 1, marks instead a second order QPT, accompanied 
by the vanishing of the energy gap in the thermodynamic 
limit. Here the fidelity vanishes only in the thermody- 
namic limit |16(. 

Finally, we observe a similar drop in the fidelity also 
at 7 = 0, for any value of pi. Here the model is exactly 
solvable and, as mentioned above, the energy gap in the 
thermodynamic limit is given by fi — 1. In spite of the 
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gapfulness of the spectrum, we identify the fidelity drop 
with a higher order QPT. Indeed, the line 7 = shows 
the same finite-size scaling properties as the transition 
far/*=l, M > 1. 




FIG. 1: GS fidelity in the complete graph for L = 400. The 
plotted function is jF m ; n = min[J-(Z, Zg^), J~(Z, Zs-y)] (see 
text), with S/j, — ^7 = 0.1. Note the symmetry with respect 
to the 7 = axis. 

Conclusions- In this paper we analyzed general sys- 
tems of free-fermions. We have showed that the informa- 
tion about the zero-temperature phase diagram of these 
linear system is encoded in the polar decomposition of 
the L x L matrix of their coupling constants; the ana- 



lyticity properties of this decomposition dictate those of 
the ground state. We have given an explict expression for 
the fidelity between different ground states in terms of the 
unitary parts of the corresponding polar decompositions 
and analyzed the ground state phase diagram in terms 
of it. In particular we showed a connection between the 
gaplessness of the many-body energy spectrum and the 
fidelity drops. Finally, to exemplify the general results 
of the paper, we have presented a preliminary study of 




-2-1012 -2-1012 

FIG. 2: Phase diagram in the fi-y plane for L even. Solid and 
dashed lines correspond to first and higher order QPTs respec- 
tively. Left panel: thermodynamic limit L — » oo. Right panel: 
phase boundaries of the first order QPT for L = 2, 4, 6, oo. 

the ground-state diagram of a free-fermi system over a 
totally connected graph. What appears to be rather re- 
markable in this fidelity approach lies in its universality 
and purely (Hilbert space) geometrical nature : no a pri- 
ori understanding of the nature of the quantum phases 
separated by critical points is required. 
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